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We report the observation of the A g -symmetric one-magnon Raman peak in the magnon Bose- 
Einstein condensation phase of TICUCI3. Its Raman shift traces the one-magnon energy at the 
magnetic V point, and its intensity is proportional to the squared transverse magnetization. The 
appearance of the one-magnon Raman scattering originates from the exchange magnon Raman 
process and reflects the change of the magnetic-state symmetry. Using the bond-operator represen- 
tation, we theoretically clarify the Raman selection rules, being consistent with the experimental 
results. 
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I. INTRODUCTION 



Currently, many physicists are examining the Bose- 
Einstein condensation (BEC) of atoms in ultracooled 
dilute gases, and in particular, the BEC of magnons. 
The latter, which is the magnetic-field induced quan- 
tum phase transition to the magnon BEC phase, has 
been reported in S = 1/2 antiferromagnets with a spin 
gap, such as KCuCl 3 , TlCuCla^M BaCuSi 2 6 £ and 
Pb2V30t). 6 The change of the magnon dispersion rela- 
tion in TICUCI3 through the magnon BEC phase tran- 
sition at H c ~ 6 T has been observed by inelastic neu- 
tron scattering^ and has been explained using the bond- 
operator representation^ One of the characteristic fea- 
tures of the magnon BEC phase is the formation of mass- 
less excitation, i.e., the Goldstone mode at the magnetic 
r point, indicating the spontaneous breaking of the con- 
tinuous symmetry. However the details of the magnon ex- 
citations, especially their symmetries, have not yet been 
established. Raman scattering is a powerful tool to study 
phase transitions. Because the magnon Raman process 
is sensitive to the symmetries of the ground and ex- 
cited states j£ Raman-scattering measurement above H c 
presents great potential to study the change of the ground 
and excited states through the magnon BEC phase tran- 
sition. 

This paper reports the observation of one-magnon Ra- 
man scattering originating from changes of the ground 
and excited states through the magnon BEC phase tran- 
sition. This study focused on TICUCI3 where the magnon 
excitations and magnetic parameters below and above H c 
have been studied in detailj££ First, we show our exper- 
imental results above H c . We then construct the micro- 
scopic theory of one-magnon Raman scattering in the ex- 
change magnon Raman process using the bond-operator 
representation, which can explain the experimental re- 
sults clearly. Based on our results, the Raman selection 



rule will be clarified. 



II. EXPERIMENTS 

Single crystals of TICUCI3 were prepared by the verti- 
cal Bridgman method^ The 5145-A line of Ar + -ion laser 
polarized along the (201) axis was incident on the (010) 
cleavage surface. We set the samples in the cryostat un- 
der the dried N2 or He gas atmosphere, because the sam- 
ple was easily damaged by moisture in air. We placed 
the microscope in the vacuum chamber of superconduct- 
ing magnet in order to collect the scattered light effec- 
tively. This enabled us to select good surface positions 
of crystals and the effects of the direct scattering in the 
low-energy region were avoided. Magnetic fields of up to 
10 T were applied nearly parallel to the (010) axis. The 
effect of a weak component of magnetic field along the 
(201) axis due to the experimental setting is negligible 
because the effect of the anisotropic g tensor along these 
directions is smallJ^ 



III. RESULTS 

Figure [1] compares the low-temperature Raman spec- 
trum at 9 T (above H c of TlCuCl 3 ) to that at T. At 
T, we observed several sharp phonon peaks superim- 
posed on the two-magnon Raman band extending from 
11 cm -1 , i.e., twice the energy of the magnetic gap^ to 
about 150 cm -1 . Thi s sp ectrum is consistent with the 
results of refs. M and Qj. At 9 T, we observed the new 
Raman peak with a Lorentzian lincshapc, called PI, at 
20 cm -1 . No other significant change was observed. PI 
excitation has the A g symmetry, which is obtained with 
the following procedures. For the incident laser polarized 
along the (201) axis, -Ein/ / (201 ) , we measured the scat- 
tered light with polarization E sc which is rotated from 
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FIG. 1: Polarization characteristics of Raman spectra at 9 
T. The inset shows the integrated Raman intensity of the 
phonon peak at 137 cm -1 . The arrows indicates the Raman 
peaks coming from the B g phonons. The unpolarized Raman 
spectrum at T is also shown. 



-Ein with an angle 8. The Raman intensities are normal- 
ized so that the 137-cm" 1 Ag-symmetric phonon peaks 
in each spectrum, for which the ^-dependence is shown in 
the inset of Fig. [T] have the same intensity. The Raman 
intensity from the quasiparticles with A g symmetry, in- 
cluding PI, is ^-independent in this plot while those with 
B g symmetry indicated by arrows increased with increas- 
ing 8, as shown in Fig. Q] 

Figure [2] shows the detailed magnetic-field dependence 
of Raman spectra at 1.9 K. At T, the 25- and 32- 
cm _1 phonon Raman peaks, called P2 and P3, respec- 
tively, are superimposed on the two-magnon Raman band 
starting at 11 cm -1 . P2 and P3 have the A g symmetry 
as well as PI, as shown in Fig. [T] The Raman spec- 
tra below 5 T are magnetic-field independent. Above 
7 T, we clearly observed that the frequency and inten- 
sity of PI strongly depended on the applied magnetic 
field, as denoted by the hatched areas in Fig. [2j of which 
the details will be explained later. Around 6 T, the in- 
crease of the Rayleigh scattering around cm -1 suggests 
the quasielastic (or critical) light scattering reflecting the 
large magnetic specific heat around H c as observed in 
several antiferromagnets or spin-Peierls system i 13 i 14 i 15 To 
discuss the quasielastic light scattering quantitatively, 
Raman-scattering measurements in the anti-Stokes re- 
gion are necessary. 

The lineshapes of PI, P2, and P3 are well described 
by three Lorentzian curves superimposed on the back- 
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Raman Shift (cm 1 ) 

FIG. 2: Magnetic-field dependence of Raman spectra in 
TICUCI3 at 1.9 K. The fitting curves are superimposed on 
the experimental data above 7 T. The details of the fitting 
curve (solid curves) together with the background generated 
by the two-magnon Raman band (dashed lines) are given in 
the text. The hatched areas show the component of the Ra- 
man intensity generated by PI. 

ground: 

t, \ (n + \)kfujTi 

7 M = 2^ r~2 2 \ 2 _i_ 1 f \2 + background , (1) 

where hi, Hu>i, and Ti indicate the Raman coupling coef- 
ficient, the energy, and the halfwidth of Pi, respectively. 
Here, the Bose factor (n + 1) can be treated as unity 
because the temperature is much lower than the ener- 
gies of quasiparticles. The background generated by the 
two-magnon Raman band peaking around 50 cm -1 was 
assumed to be a linear function, as shown by the dashed 
lines in Fig. [2] The calculated curves reproduced the 
observed data well. We show the Raman intensity gen- 
erated by PI (the term related to the subscript i = 1 in 
eq. ([T])) on the linear background as the hatched area in 
Fig. O Around H c , we could not distinguish PI from the 
two-magnon Raman band because of its weak intensity. 
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Magnetic Field (T) 

FIG. 3: (a) Magnetic-field dependence of one-magnon (filled 
circles) and phonon (empty circles) energies together with 
that of one-magnon energies calculated in refs. |8| and 

HE (b) 

Magnetic-field dependence of squared Raman coupling coef- 
ficient (filled circles). Squared transverse magnetization M^ y 
(empty squares, ref. HI) and squared longitudinal magneti- 
zation Ad? (filled squares, ref. uM ) are also plotted with the 
calculated values for them (a solid and a dashed curves for 
Ml y and M%, respectively) £ 



Figure 02a) shows the peak energies fuox and hu>2 
as functions of magnetic field together with the calcu- 
lated one-magnon energy E ga (Q) (a — — , 0, +) with 
the wavevector Q = (0, 0, 27r)j&i& where the magnetic 
gap is closedi^iiLi One can see that hio\ below 10 T 
agrees with E g —(Q) within experimental accuracy. Fig- 
ure EHb) shows the squared Raman coupling coefficient 
k 2 = kf/k 2 , which is proportional to the integrated Ra- 
man intensity of PI. Here, we normalized k 2 , which is 
proportional to the area hatched in Fig. 03 by k 2 to cor- 
rect errors due to the small deviations of optic alignment. 
The errors of k 2 are similar to the symbol size in Fig. 
OUJb). For comparison, we show the magnetic- field de- 
pendences of squared transverse magnetization M 2 y (ref. 
H) and squared longitudinal magnetization M 2 (ref. fill ) 
together with their calculated values 4 One can see that 
the magnetic-field dependence of Raman intensity is well 
scaled to the former. 



IV. DISCUSSION 

First, we discuss the origin of PI. Judging from the 
polarization characteristics, the scattering process of PI 
comes from the exchange magnon Raman scattering, as 
proposed by Fleury and Loudon^ which usually cre- 



ates the broad two-magnon Raman band with A g sym- 
metry reported at zero magnetic field i 11 ' 12 Because the 
two-magnon Raman band at 9 T is almost the same 
as that at T, as seen in Fig. [TJ we do not need to 
consider the increase of the magnon-magnon interaction 
which may cause the two-magnon Raman band with a 
nearly Lorentzian lineshape— ^ as well as the forma- 
tion of the two-magnon bound stated! PI does not orig- 
inate from a three-magnon process, where the thermally 
excited triplets play an essential role i 22 i 23 The three- 
magnon Raman intensity, which is proportional to n at 
low temperatures, should be negligibly small at 1.9 K. 
Because PI has a Lorentzian lineshape, the one-magnon 
Raman scattering can be considered as the origin of PI. 
Hereafter, we consider the detail of the magnon Raman 
process using the bond-operator representation and clar- 
ify that the one-magnon Raman scattering from the ex- 
change magnon Raman process becomes possible in the 
magnon BEC phase. 

The effective Raman operator TZ in the above- 
mentioned exchange magnon Raman process has a form 
given by the isotropic Heisenberg-type exchange interac- 
tion between spins Si and Sj: 

n = J2 U ^ = F h] {E in ■ r l3 ){E sc ■ r ZJ )S t ■ Sj , (2) 

where r.y indicates the position vector between Si and 
Sj and the sum runs over all the interacting spin pairs. 
Here, r — r/\r\, E = E/\E\, and the coefficient F^j 
depends on the pathway of interaction between Si and 
Sj. TZ depends on the experimental setting through the 
(E^rij ) (i£ sc -r y ) term. The matrix element of TZ between 
the initial state |i) and the final one |f) is called Raman 
tensor. The magnon Raman intensity is given as 

I (in ' sc) M« |£ in | 2 ^|(i|ft|f)| 2 ^- Wif ) , (3) 
If) 

where fnon is the excitation energy between the states |i) 
and |f). 

We note here that TZ is always written using the 
pure singlet operator s k and the triplet operators tka 
{a = — , 0, +), which annihilate the triplets with S z = a. 
We need to rewrite TZ using the creation and annihila- 
tion operators of eigenstates in the magnon BEC phase. 
As discussed by Matsumoto et al.j^& the following trans- 
formed operators based on the bond-operator representa- 
tion characterize the magnon excitations in the magnon 
BEC phase: 

s k =ua k -vb k +Q+, 

tk+=vfa k -Q+ufb k+ -gb k -, 

t k o~b k0l 

t k -=vga k - Q +ugb k+ + fb k -, 

where the momentum-independent real-number parame- 
ters it, v, /, and g satisfy f 2 +g 2 = u 2 +v 2 = 1. Below H c , 
v = and one can obtain the simple relations a k = s k 



and bk a — tk a . Above H c , v ^ and the operators ak-Q 
and b^± are linearly combined. The mixing of bf.± and 
b^_ k± is treated using the Bogoliubov transformation, as 
will be shown in detail later. It should be noted that the 
ground state does not include the bko state, indicating 
that the E g0 (Q) mode is not Raman-active. 

We consider the Raman operator lZ d associated with 
the intradimcr interaction, which can be written in the 
reciprocal lattice space as 

K d = F d (E in -d)(E sc -d)J2[-ls{s k + X~*L**«) 

k \ a ) 



F d (E in -d){E sc -d) 



— — U a 



-uva(bQ-\ 



+ 



E 

k 



-fi-bk- 



\ b lo h kO 



(5) 

where d indicates the position vector between two spins 
forming a dimcr, which is almost parallel to the (201) 
direction, and F d originates from Fij. Here, we used 
the facts that the operator ak can be treated as a uni- 
formly condensed mean-field parameter aSk,o 

(ref. M) 

and b^_Q + = 6q+ at the magnetic V point. The term 

uva(bQ + + 6q + ) in eq. ([5]) gives the momentum selection 
rule of one-magnon Raman scattering. The parameter uv 
indicates the appearance of one-magnon Raman scatter- 
ing only above H c and the Raman intensity proportional 
to Ml y (see eq. (5) of ref. which is consistent with 
the observation, as shown in Fig. OUb). The terms a 2 
and b\, a bk a in eq. do not give one-magnon Raman 
scattering. 

The quadratic terms of the magnetic Hamiltonian H± 
in ref. [3 can be diagonalized using the a k bosonic op- 
erators which annihilate the E g ±(k) modes, respectively. 
The ground state in the magnon BEC phase is the vac- 
uum state for the a k operators. Using the bosonic com- 
mutation relations of a k , we obtain the inverse Bogoli- 
ubov transformation as 



f bk -) 
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bk+ 










\ 



u 



+* 
-fe- 



— * -4-sft 

u -k+ u -k- 

— * 

-V 1 ; 



J -k~ 

-k+ 



-k- 



l k+ 



L k- 
l k+ , 



(6) 



The one-magnon term in eq. ([5]) can be rewritten as 
b Q+ + b Q+ = ( U Q+ ~ v Q+T a Q + ( u 



l Q+ 



(7) 

indicating that both the E g+ (Q) and E g -(Q) modes are 
symmetry-allowed. These modes can be described as the 
mixed amplitude and phase modes which are related to 
the spin correlation functions along the x and y directions 
in Fig. |4l respectively. The amplitude mode changes 
the amplitudes of M xy without changing their directions 
whereas the phase mode is described as uniform rota- 
tions of M xy . Both of these modes have A g symmetry, 
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FIG. 4: Changes of spin polarizations for the amplitude and 
phase modes, corresponding to the ionic-vibration patterns 
in phonons. The solid and dashed lines indicate the dimers 
at the corner and center of the chemical unit cell denoted by 
parallelograms, respectively. Only the Cu 2+ sites are shown 
with the symbols of inversion centers and screw axes. The 
direction of the transverse magnetization x in the magnon 
BEC phase is shown by a bold arrow with M xy £ The direction 
y is perpendicular to x and b. 



i.e., these are intrinsically Raman-active because the con- 
tinuous rotational symmetry is broken above H c . This is 
one of the most distinguishing characteristics of the mag- 
netic excitations in the magnon BEC phase. In systems 
of density waves^i the Goldstone mode (phase mode), 
which corresponds to the continuous translational oper- 
ation, is IR-active. 

Let us show that the phase mode does not give the fi- 
nite Raman intensity. When the phase mode is the Gold- 
stone mode, i.e., for E g+ (Q) = 0, TL± can be diagonalized 
by 













A Q± 




-Aq± 


e Q- 


\Aq± 




Aq- 




= 0, 



(8) 



-eg-/ \v+J 
where the definitions of eq^ and Aq^ (/? = - 



, ±) are 
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given in ref. H|. We can reduce eq. ([8]) to the form of 

fe Q+ + A Q+ e Q± + A Q± Y u+ + - v+\ _ 

\e Q± + Aq± e Q _ + A Q -){u+_ v%_) U " {J) 

Because the matrix in eq. ^ is invertible, we find that 
Uq + — v~q + = 0, i.e., the Raman intensity for the phase 
mode is zero although this mode is symmetry-allowed. 
This result indicates that only the spin correlation func- 
tion along x is detectable with a factor M 2 y by the first- 
order Raman scattering. 

In case of TICUCI3, the anisotropic exchange interac- 
tion gives a small magnetic gap E g+ (Q) rj 1.7 cm -1 ^ 
In this case, the E g+ (Q) mode is the mixed amplitude 
and phase modes and it gives a finite one-magnon Ra- 
man intensity. In our measurements, however, we could 
not detect it because of the strong direct scattering at 
cm -1 . It is worthwhile to consider the Raman scat- 
tering from the E g+ (Q) because this mode is thermally 
populated at 1.9 K. The transition to the ground state 
is the one-magnon anti-Stokes Raman scattering, which 
cannot be detected in our measurements, as stated above. 
The transition to the E g _(Q) mode is obtained from the 

terms ot~^ in Raman tensor. Substituting and 

in eq. ^ for bk± and w k± in eq. (JSJl, IZd contains the 

terms 




(10) 

This result indicates that the transition from the E g+ {Q) 
state to the E g _ (Q) one may be detected as a part of the 
two-magnon Raman band and its intensity is expected to 
have no drastic change, at least below 10 T, because v 2 
below 10 T is very small4^ Actually, the profile and in- 
tensity of the two-magnon Raman band at 9 T are almost 
similar to those at T, as shown in Fig. [TJ 

When we consider one-magnon Raman scattering 
caused by the interdimer interaction, we can substitute 
the expectation value for one of the spin operators in eq. 
©: 

Tl = F itj (E ia • rl 3 )(E sc • r^Si ■ (Sj) . (11) 

The Raman tensor from the interdimer interaction also 
contains the terms (6q+ + &q + ), suggesting that the one- 
magnon Raman scattering from the E g _ (Q) mode can be 
detected above H c . Because (S) — M xy , one can expect 
that the Raman intensity is also proportional to M 2 y . 
The precise analytic form of the Raman tensor coming 



from the interdimer interaction has been established in 
our recent letter in case of the pressure-induced magnon 
BEC phase at zero magnetic field.™ The magnetic field- 
induced magnon BEC case will be published elsewhere^ 

We point out that our theory for the appearance 
of the one-magnon Raman scattering is applicable to 
the pressure-induced magnon BEC phase transition in 
TlCuCi3i & 27 ' 28 In case of pressure- induced magnon BEC 
phase transition, the pure amplitude mode, which is the 
longitudinal spin-wave mode coupled only with the spin 
correlation function along the x direction in Fig. 31 is 
expected to be observedi 25 ' 29 

The two-magnon Raman band is also interesting as 
well as the appearance of the one-magnon Raman peak 
which is the main purpose of this paper. In case of two- 
magnon Raman scattering, both IZd and the Raman ten- 
sors generated from the interdimer interactions, which 
create the magnon pair with the zero total momentum, 
play an essential role. At present, it was difficult to cal- 
culate the lineshape of the two-magnon Raman band be- 
cause the values of Fij, which the two-magnon Raman 
spectrum is sensitive to, cannot be obtained directly. 



V. CONCLUSION 

We have assigned the origin of the Raman peak appear- 
ing above H c in TICUCI3 to one-magnon Raman scat- 
tering, which comes from the exchange magnon Raman 
process. This is based on (1) the Lorentzian lineshape 
of the peak, (2) its Raman shift tracing E g -(Q), (3) its 
polarization characteristics, i.e. this one-magnon Raman 
scattering is j4 g -symmetric as well as the second-order 
magnetic Raman scattering, and (4) the observation that 
the peak's Raman intensity is proportional to M 2 y . Us- 
ing the bond-operator representation, we calculated the 
Raman intensity to clarify the Raman selection rule of 
one-magnon Raman scattering in the exchange magnon 
Raman process. The intensity of the one-magnon Ra- 
man scattering is related to the spin correlation function 
along the direction of M xy , i.e., the x direction in Fig. 
4. And therefore the E g ±(Q) modes with the finite exci- 
tation energies are A g -symmetric and Raman-active. In 
the isotropic limit, the Goldstone mode for E g+ (Q) = 0, 
which is related to the spin correlation function along the 
y direction, is Ag-symmetric but has no Raman inten- 
sity. The E g o(Q) mode and the magnetic excitation at 
the chemical T point are Raman-inactive. The change of 
the ground and excited states through the magnon BEC 
phase transition can be detected via the appearance of a 
new one-magnon Raman peak from the E g _ (Q) mode in 
the magnon BEC phase. 



* Electronic address: kuroe@sophia.ac.jp 
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